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Quantum  repeaters  create  long-distance  entanglement  between  quantum  systems  while  overcoming  difficul¬ 
ties  such  as  the  attenuation  of  single  photons  in  a  fiber.  Recently,  an  implementation  of  a  repeater  protocol 
based  on  single  qubits  in  atomic  ensembles  and  linear  optics  has  been  proposed  [Duan  et  ai,  Nature  (London) 

414,  413  (2001)].  Motivated  by  rapid  experimental  progress  towards  implementing  that  protocol,  here  we 
develop  a  more  efficient  scheme  compatible  with  active  purification  of  arbitrary  errors.  Using  similar  resources 
as  the  earlier  protocol,  our  approach  intrinsically  purifies  leakage  out  of  the  logical  subspace  and  all  errors 
within  the  logical  subspace,  leading  to  greatly  improved  performance  in  the  presence  of  experimental  ineffi¬ 
ciencies.  Our  analysis  indicates  that  our  scheme  could  generate  approximately  one  pair  per  3  min  over 
1280  km  distance  with  fidelity  (F3S78%)  sufficient  to  violate  Bell's  inequality. 

DOI:  10. 1103/PhysRevA.76.0 12301  PACS  number(s):  03.67.Hk,  42.50. -p,  42.50.Gy 


I.  INTRODUCTION 

Quantum  communication  holds  promise  for  the  secret 
transfer  of  classical  messages  as  well  as  forming  an  essential 
element  of  quantum  networks,  allowing  for  teleportation  of 
arbitrary  quantum  states  and  violations  of  Bell’s  inequalities 
over  long  distances  [1],  While  experimental  and  even  com¬ 
mercial  implementation  of  simple  quantum  communication 
protocols  are  well-established  [2,3],  extending  these  tech¬ 
niques  to  distances  much  longer  than  the  attenuation  length 
of  optical  fiber  remains  a  challenging  goal  due  to  exponential 
attenuation  of  transmitted  signals.  Quantum  repeaters  [4-6] 
overcome  the  exponential  time  overhead  associated  with  fi¬ 
ber  attenuation  and  other  errors  by  using  a  quantum  memory 
and  local  quantum  computation. 

Several  promising  avenues  for  quantum  repeater  imple¬ 
mentation  include  both  atomic  ensembles  [7]  and  using  few 
qubit  quantum  computers,  such  as  neutral  atoms  in  cavity 
QED  [8,9],  ion  traps  [10],  and  solid-state  single  photon  emit¬ 
ters  [11].  Experimental  progress  [12-14]  towards  realization 
of  the  Duan-Lukin-Cirac-Zoller  (DLCZ)  protocol  [7]  has 
been  especially  rapid,  with  many  building  blocks  demon¬ 
strated  in  the  laboratory.  The  experimental  challenge  is  now 
shifting  towards  the  realization  of  scalable  quantum  repeater 
systems  which  could  yield  a  reasonable  communication  rate 
at  continental  distances  (a  1000  km).  Thus  the  DLCZ  proto¬ 
col  should  be  examined  and  adapted  to  practical  experimen¬ 
tal  considerations,  allowing  one  to  remove  imperfections 
such  as  the  finite  efficiency  of  retrieval  and  single-photon 
detection  and  fiber  length  fluctuations.  Our  approach  extends 
the  DLCZ  protocol,  keeping  the  experimental  simplicity  of 
the  original  scheme  while  avoiding  fundamental  difficulties 
due  to  these  expected  experimental  imperfections. 

This  paper  is  organized  as  follows.  In  Sec.  II  we  will 
review  the  DLCZ  protocol  and  describe  our  approach  which 
uses  a  basis  to  encode  each  qubit.  Section  III  compares  both 
the  DLCZ  protocol  and  our  approach  in  the  presence  of  im¬ 
perfections.  Section  IV  estimates  the  time  scaling  of  our  ap¬ 
proach  and  compares  three  specific  implementations.  Section 
V  summarizes  our  results. 


H.  ATOMIC-ENSEMBLE-BASED  QUANTUM  REPEATERS 
A.  DLCZ  protocol:  Review 

The  DLCZ  protocol  [7]  starts  with  entanglement  genera¬ 
tion  (ENG)  by  counting  the  interfering  Stokes  photons  scat¬ 
tered  from  a  pair  of  distant  atomic  cells  x  and  y.  This  gener¬ 
ates  an  entangled  state 

\Qx,y=(Sl  +  e^Sl)/^2\vaC)x,y,  (1) 

with  .S']  and  y  the  creation  operators  of  spin-wave  modes  in 
the  two  cells,  respectively,  and  0  the  phase  difference  be¬ 
tween  left  and  right  channels  for  Stokes  photons  [15].  Then, 
entanglement  connection  (ENC)  is  performed  on  two  pairs  of 
entangled  atomic  cells  |£01)^>y  and  \^)XoyR,  obtaining  a 
further  separated  entangled  pair  \£,li1]+lf,^}X]  yR  probabilistically. 
The  ENC  step  provides  built-in  purification  against  many 
imperfections — photon  loss,  atomic  excitation  loss,  and  dark 
counts.  In  the  final  step,  postselection  is  used  to  obtain  an 
effectively  polarization  entangled  state 

|'P™£>  =  e^(Slsl2  +  Slfy/^2  |vac>  (2) 

from  two  parallel  pairs  |  ^xvy^x2,y2,  which  overcome 
static  phase  errors  (time  independent  0’s). 

There  are  two  important  merits  of  the  DLCZ  protocol. 
First,  it  has  intrinsic  purification  of  errors  due  to  photon  loss 
(in  the  fiber,  the  quantum  memory,  and  the  photon  detector) 
and  significantly  relaxes  the  experimental  requirement  for 
quantum  repeater.  In  addition,  the  time  scaling  of  the  DLCZ 
protocol  is  always  subexponential  and  very  close  to  polyno¬ 
mial  when  the  retrieval  and  detection  efficiency  is  high. 
However,  the  DLCZ  protocol  does  not  purify  all  kinds  of 
errors.  For  example,  time  dependent  0’s  (due  to  fiber  length 
fluctuation)  induce  phase  error,  which  cannot  be  taken  out  as 
a  common  factor  in  Eq.  (2),  since  the  two  pairs  of  entangled 
atomic  cells  are  not  produced  at  the  same  time.  Such  a  phase 
error  is  accumulated  and  doubled  after  each  level  of  ENC.  In 
addition,  combined  photon  loss  during  ENG  and  ENC  may 
also  induce  a  phase  error  not  purified  by  the  DLCZ  protocol. 
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FIG.  1.  (Color  online)  Repeater  components:  (a)  entanglement  generation  (ENG)  and  (b)  entanglement  connection  (ENC);  indicated 
operations:  retrieve  bc  and  ac  (additional  45°  rotations  only  for  the  first  level),  join  on  polarizing  beam  splitter  (PBS),  detect  in  ±  basis 
conditioned  on  one  photon  per  output,  and  finally  adjust  the  phase,  (c)  Entanglement  purification  (ENP);  indicated  operations:  retrieve  cq.Zq 
and  a2,b2  (additional  45°  rotations  to  purify  phase  error),  interfere  alta2  on  PBS  (same  with  bt  ,b2),  restore  a3,i3  conditioned  on  single 
photon  at  a4  and  b4  respectively,  and  finally  adjust  the  phase. 


Furthermore,  the  DLCZ  protocol  (dashed  line  in  Fig.  2)  still 
has  a  significant  time  overhead  for  long  distances  because  of 
the  superpolynomial  scaling  in  the  presence  of  realistic  im¬ 
perfections.  For  instance,  nonideal  retrieval  and  detection  ef¬ 
ficiency  (rj=  7} retrieve  X  ^detection  <  1)  during  ENC  introduces 
a  large  vacuum  component,  suppresses  the  success  probabil¬ 
ity  of  later  ENC,  and  consequently  slows  down  the  protocol. 
Motivated  by  these  issues,  we  will  extend  the  DLCZ  proto¬ 
col,  mitigating  the  above  errors. 

B.  Different  approach 

We  now  consider  a  different  approach  in  which  two 
atomic  cells  are  used  at  each  node  a,  labeled  ( a,H )  and 
( a,V ),  to  store  one  qubit,  a.  The  qubit  is  defined  as  one 
single  spin- wave  excitation  shared  between  two  cells: 

{|tf>fl  =  <ff|vac>,|V>n  =  <y|vac>}.  (3) 

When  the  stored  spin  waves  are  converted  back  into  photons, 
the  photons  have  a  polarization  ( H  or  V)  consistent  with  that 
stored  in  the  originating  cell.  This  qubit  basis  allows  projec¬ 
tive  measureme_nts  along  any  qubit  states,  e.g.,  |±) 

=  {\H)a±\V)a)l \2,  using  linear  optical  operations  and  photon 
counting  [16].  We  will  show  that  in  this  logical  basis  it  is 
possible  to  perform  entanglement  purification  (ENP)  [6]  to 
reduce  errors  within  the  logical  subspace,  including  phase 
fluctuation.  Since  ENP  can  suppress  errors  within  the  logical 
subspace  which  occur  with  probability  q  to  0{q2),  only  a  few 
ENP  levels  are  needed  to  obtain  high  fidelity  entanglement. 

We  now  describe  our  procedures  for  ENG,  ENC,  and 
ENP.  ENG  [Fig.  1(a)]  is  similar  to  that  of  the  DLCZ  proto¬ 
col,  but  here  two  parallel  entangled  pairs  are  generated  be¬ 
tween  a  and  b: 


l^£A,G)a,i  -  \£<t)(a,H)(b,H)\£cf^(a,V)(b,V) 

=  e^(\H)a\V)b  +  |  V)a\H)b)  +  |ffV»ac>* 

+  e2'^|vac)a|//V%.  (4) 

The  entangled  states  are  prepared  in  the  quantum  memory,  so 
no  simultaneity  is  required  for  creating  the  two  states  com¬ 
prising  lip™0).  For  small  excitation  probability  pc,  the  whole 
generation  only  takes  time  0(l/pc),  in  contrast  to  0(l/pf) 
for  schemes  requiring  simultaneity,  e.g.,  coupling  between  a 
trapped  atom  and  photon  [10]  or  parametric  down  conver¬ 
sion  [17].  Errors  from  multiphoton  events  occur  only  with 
probability  /?2,  and  are  considered  in  a  later  analysis  of  im¬ 
perfections. 

The  first  level  of  ENC  converts  two  |'P£3VG)  states  (one 
between  aL  and  bc,  the  other  between  ac  and  bR)  into  polar¬ 
ization  entangled  states  \^)ab=\H)a\lH)bR+\V)a^V)bR.  Only 
4  out  of  the  16  terms  in  the  Schmidt  decomposition  of 
\^rENG)a[b(\1^rENG)acbR  have  any  contribution  to  the  output 
state;  the  remainder  are  eliminated  by  projective  measure¬ 
ment  during  ENC,  reducing  the  probability  of  success  for 
ENC  from  1/2  to  1/8.  At  higher  levels  of  ENC,  the  opera¬ 
tions  correspond  to  standard  entanglement  swapping  [4,5], 
where 

1^*  ) aLbC  ®  ^  ) aCbR  ) aLbR 

leads  to  an  entangled  pair  between  L  and  R  with  probability 
1/2,  as  detailed  below. 

The  procedure  for  ENC  is  illustrated  in  Fig.  1(b).  First, 
the  spin  waves  stored  in  qubit  bc  and  ac  are  retrieved  into 
photons.  At  the  lowest  level  of  ENC,  the  polarization  of  the 
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TABLE  I.  Entanglement  connection  procedure  applied  to  I®*) 
inputs  for  entangled  pairs  between  aL  and  bc,  and  ac  and  bR.  For 
clarity,  we  introduce  |0)  =  |H)  and  1 1 )  =  |  V)  to  represent  logical 
states  (i.e.,  states  with  exactly  one  excitation),  | HV)  for  nonlogical 
states  with  two  excitations,  and  |vac)  for  states  with  no  excitation 
(sometimes  omitted).  We  assume  H  photons  pass  through  and  V 
photons  are  reflected  at  the  middle  PBS. 


Operation 

Transform  of  |<P±>az6c|<l>±>aci,R 

Retrieve 

bc,ac 

|0000)±  |001 1)±  1 1 100>  +  |  1 11 1) 

Transform 
bc,ac  at  PBS 

|0000)+|1111) 

±|0>at(|tfV>flc)|l>6K 

±\l)aL(\HV>bc)\0)bR 

One  photon 
per  mode 
(p= 0.5) 

|0000)+|1111) 

Detect  in  ±, 
results  m,m’ 

W)aLbR+mm'\n)aipR 

Phase  shift 

mm1 

1 

photons  is  rotated  45°.  The  rotations  transform  \HV)„c(or  bc> 
into  (\HH)-\ VV))ac(OI bc),  because  for  bosonic  fields  [18] 

44  45  (4 + 4)  (4 - 4)  =  (44 - 44)  ■ 

Thus  after  the  polarizing  beam  splitter  (PBS)  there  will  be  at 
least  two  photons  at  one  output.  For  incoming  state 
\'&ENG)aLbc<&\'irENG)acbR,  all  seven  terms  containing  two  ex¬ 
citations  in  at  least  one  pair  of  cells  in  the  center  repeater 
node  (such  as  |#V)flc(or6c))  do  not  contribute  to  the  click 
patterns  with  one  photon  at  each  output.  Five  terms  contain¬ 
ing  two  excitations  in  one  of  the  left  or  right  repeater  nodes 
[e.g.,  | HV>  (od>  j]  have  at  most  one  excitation  retrieved  from 
bc  and  ac,  which  is  insufficient  to  give  two  clicks.  Therefore 
only  the  four  terms  remaining  can  give  the  correct  photon 
detector  click  patterns. 

For  all  levels  of  ENC,  the  photons  are  then  joined  on  the 
middle  PBS  and  the  number  of  photons  at  two  outputs  are 
counted  in  the  {|+),|— )}  basis.  With  probability  50%,  there  is 
one  photon  at  each  output,  and  the  connection  is  successful; 
otherwise  the  process  is  repeated.  If  the  two  photons  have 
orthogonal  polarizations,  a  bit  flip  a\H)  +  f3\V)—>  a|V> 
+  /3\ H)  is  applied  to  aL  [19].  At  higher  levels  of  ENC,  where 
the  45°  rotations  are  not  necessary,  the  bit  flip  is  replaced 
by  the  phase  flip  a\H)  + /3\V) ^  a\H) - (3\V) ,  as  detailed  in 
Table  I  [20]. 

The  third  component  is  ENP  [Fig.  1(c)]  which  obtains  a 
high  fidelity  entangled  pair  from  two  pairs.  Our  procedure 
uses  polarization  entangled  photons  and  is  similar  to  recent 
experimental  investigations  [17].  During  entanglement  puri- 


TABLE  II.  Bit  entanglement  purification  acting  on  |<F+)|1P+). 
The  same  notation  as  Table  I  is  used.  For  such  given  input  states, 
the  click  pattern  never  matches  the  right  pattern,  and  bit-error  of  a 
single  input  qubit  can  be  filtered  completely.  The  remaining  steps  to 
preserve  the  desired  Bell  states  are  detection  of  a4,b4  in  ±  basis,  a 
phase  shift  n  n’  based  on  the  measurement  results,  and  storing  of 
the  traveling  photons  (a2,b3)  into  atomic  ensembles. 


Operation 

Transform  of  \$>n)aib^n  )a^ 

Retrieve 

ctlbl,a2b2 

|0001)+;i|  1 101)  +  rc'|0010)  +  ;m'|ll  10) 

Interfere  a\,a2 
on  PBS;  same 
with  bl,b2 

|00)a3a4(|//V%4) 

+n(|tfV>a3)|ll)i3i4 

+n'(\HV)aA)\00)b}bA 

+nn'\U)a^(\HV)b) 

One  photon 
per  lower 
output 

0 

fication  of  bit  errors  (bit-ENP),  the  qubits  from  two  parallel 
pairs  pc,l  bl  and  plh  h^  are  retrieved  from  the  quantum  memory 
and  joined  at  PBSs.  The  photons  for  two  upper  outputs  are 
stored  into  quantum  memory  a3  and  7>3.  The  photons  for  the 
lower  outputs  a4  and  b4  are  counted  in  {| +),|  — )}  basis.  With 
probability  50%,  there  is  exactly  one  photon  at  each  lower 
output,  and  the  purification  is  successful;  otherwise  two  new 
pairs  are  created  by  restarting  the  process.  If  the  two  photons 
have  orthogonal  polarizations,  a  phase  flip  is  applied  to  a3. 
An  example  of  purification  of  bit-error  is  presented  in  Table 
II.  During  purification  of  phase  errors  (phase-ENP),  addi¬ 
tional  45°  rotations  are  applied  to  the  retrieved  qubits  and  the 
bit  flip  is  replaced  by  the  phase  flip.  The  addition  of  45° 
rotations  effects  the  basis  transform  |''P+),  leading  to 

purification  of  errors  of  the  other  type.  The  truth  table  of 
phase-ENP  is  listed  in  Table  III.  Bit  (or  phase)  errors  can  be 
nonlinearly  suppressed  to  the  second  order  during  bit-ENP 
(or  phase-ENP)  [21], 

The  three  components  described  above  for  the  quantum 
repeater  protocol  only  use  atomic  cells,  linear  optics,  and 
photon  number  counting.  We  remark  that  the  duration  of  the 
retrieved  anti-Stokes  pulse  can  be  made  long  (a  I  /is)  com¬ 
pared  to  the  detector  recovery  time  by  adjusting  the  intensity 
and  duration  of  the  retrieval  pulse.  This  enables  photon  num¬ 
ber  counting  of  the  anti-Stokes  pulse  [22]. 

TABLE  III.  Truth  table  for  phase-ENP.  Each  element  gives  the 
possible  output  state  after  the  purification  operation.  (“■ — ”  for  cases 
with  no  outputs.) 


'P<72^2 

<|.+ 

\J/+ 

\J /- 

<F+ 

4>+ 

— 

r+ 

— 

— 

4>- 

— 

\J /- 

\J/+ 

\J/+ 

— 

<t>+ 

— 

\]/“ 

— 

X\ t~ 

— 

4>- 
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III.  NOISE  AND  IMPERFECTIONS 


A.  Nonlogical  errors 


We  now  examine  the  performance  of  our  new  scheme  by 
considering  the  role  of  errors,  starting  with  how  imperfec¬ 
tions  due  to  inefficiency  limit  the  protocols.  Primarily,  we 
find  that  inefficiency  takes  logical  states  into  two  types  of 
nonlogical  states — those  with  too  few  excitations  ( vacuum 
type)  and  those  with  too  many  excitations  ( multiexcitation 
type).  We  represent  these  errors  by  density  matrix  7 rvac  (a 
mixed  state  with  at  most  one  excitation  between  both  pairs  of 
cells)  or  TTmuiti  (a  mixed  state  with  at  least  one  pair  of  cells 
with  more  than  one  excitation).  The  normalized  density  ma¬ 
trix  after  mth  ENC  [Fig.  1(b)]  can  be  written  as 


Jm)  _  Am)  ( m )  ,  Am)  On)  Am)  (m) 

•  a^,b ^  Jr  logic r* logic  r  vac  11  vac  r  multi '  multi  ’ 


(5) 


where  the  m -dependent  operator  p,"'j  is  the  density  matrix 
within  logical  subspace;  7 r  and  7 r" ...  also  depend  on  m: 
and  the  coefficients  p,  .  ,  p  '">,  and  p  are  the  probabilities 

r  lOglC ’  r  vaci  t'  multi  r 

for  the  logical,  vacuum,  and  multiexcitation  types,  respec¬ 


tively. 


After  the  first  level  of  ENC, 


Ji) 


L 


V  and  p(2i,rPc 


<S  1 .  We  can  demonstrate  that  these  three  probabilities  re¬ 
main  stable  (see  Appendix  B)  for  all  higher  levels  of  ENC, 
by  considering  the  unnormalized  state  after  (/w+l)th  ENC 
for  the  new  approach:  /r ™ + 11 =ph 

~(m+l)  (m+1) 

+P multi  77 multi  ’  Wlth 


— (m+ 1 )  _  — (m+ 1 )  (m+ 1 )  — (m+ 1 )  (m+ 1 ) 


logic  P logic  P vac 


77 


-im+D  _  I_(»)  Am)  r,  (m+1)  n(  n 
P  logic  2  PPlogicP  logics  +  P err, new  +  '-'\Pc)  J’ 


Am+l)  __  _ 
r  vac 


2VP^J^+0(pc)l 


^~\vP(tiJmlll+0(pc)\, 

where  the  probability  for  the  new  logical  error  from  the  mul¬ 
tiexcitation  states  (accompanied  by  photon  loss)  is 

pt:l~0-y)Pc- 

The  total  logical  error  probability  p'"^ ' 1  has  two  contribu¬ 
tions:  the  accumulated  logical  errors  from  both  input  pairs 
for  ENC,  and  the  new  logical  error 

PeJJX)  ~  2 P(err  ~  (1  “  V)Pc-  (6) 

We  can  calculate  for  the  new  scheme, 

P^~(2m-m-y)Pc-  (7) 

A  more  detailed  calculation  (see  Appendix  B),  in  which 
7 Tvac  and  7 Tmuiti  are  further  divided  into  subspaces  with  a  dif¬ 
ferent  number  of  excitations  (e.g.,  rrl:ar  is  subdivided  into 
zero-excitation  and  one-excitation  subspaces)  verifies  the  sta¬ 
bility  of  the  probability  distribution  of  p(,m\  ,  p('n\  and  pim\,- 
[Similarly,  dark  count  can  also  induce  errors  in  logical  sub¬ 
space  with  probability  ~pjark{t  -  ps),  which  is,  however, 
negligible  due  to  very  low  dark  count  probability  p^ark-i 


FIG.  2.  (Color  online)  Comparison  between  the  DLCZ  protocol 
and  the  new  scheme  (NS)  [without  active  entanglement  purification 
(ENP)].  For  each  distance,  we  optimize  over  the  choice  of  the  con¬ 
trol  parameters  (the  half  distance  between  neighboring  repeater  sta¬ 
tions,  L0,  and  the  elementary  pair  generation  probability,  pc).  With 
efficiency  p=90%  and  targeting  fidelity  F= 90%,  we  find  the  most 
efficient  implementations  to  create  the  polarization  entangled  state 
[Eq.  (2)]  for  both  the  DLCZ  (circled  black  dashed  line)  and  the  new 
scheme  (squared  blue  solid  line).  The  fiber  attenuation  length  is 
La„=  20  km,  with  no  dynamical  phase  error.  The  main  plot:  we 
show  the  relationship  between  the  (optimized)  average  creation 
time  tavg  and  the  final  distance  L  for  both  schemes,  and  the  empiri¬ 
cal  estimate  [Eq.  (11)]  of  the  time  scaling  for  the  new  scheme  (blue 
dotted  line).  Over  long  distances  (L>320km),  the  polynomial 
scaling  of  the  new  scheme  is  more  favorable  than  the  superpolyno¬ 
mial  scaling  of  the  DLCZ  protocol.  Inset:  we  plot  the  fidelities  of 
the  intermediate  distances  (Linl=  160,  320,  640,  and  1280  km)  to 
create  polarization  entangled  states  (L=1280,  F=  90%),  with  the 
optimized  choice  of  the  control  parameters  (L0,pc)  =  (80, 0.0027) 
and  (40,0.0081)  for  the  DLCZ  (tavg~  1900  s)  and  the  new  scheme 
(tavg~  380  s),  respectively.  The  optimized  choices  of  the  control 
paramters  are  detailed  in  Tables  IV  and  V. 

For  the  DLCZ  protocol,  only  two  cells  are  used  to  store 
entanglement.  Besides  the  logical  states  (single  excitation  in 
two  cells),  we  can  similarly  define  the  vacuum  states  (with 
no  excitation)  and  multiexcitation  states  (with  two  or  more 
excitations).  Contrary  to  our  approach,  the  probability  distri¬ 
bution  is  not  stable — both  vacuum  and  multiexcitation  prob¬ 
abilities  increase  with  distance  (see  Appendix  A).  The 
vacuum  probability  soon  becomes  the  dominant  term,  which 
reduces  the  success  probability  of  ENC  significantly,  result¬ 
ing  in  superpolynomial  (but  still  subexponential)  time  scal¬ 
ing  (Fig.  2).  The  logical  error  probability  for  the  DLCZ  pro¬ 
tocol  has  the  same  form  as  Eq.  (6)  up  to  the  coefficient,  but 
the  ratio  P^itJ P^iogic  (thus  p'™])  grows  with  distance  (see  Ap¬ 
pendix  A),  which  accounts  for  the  sharp  decrease  of  fidelity 
for  the  DLCZ  protocol  (see  inset  of  Fig.  2).  To  maintain  good 
final  fidelity,  the  initial  error  p^lti  (and  pc)  should  be  very 
small,  which  demands  longer  generation  time  of  an  elemen¬ 
tary  pair  for  the  DLCZ  protocol. 

In  essence,  by  requiring  at  least  one  excitation  in  the  en¬ 
semble,  our  qubit  subspace  is  automatically  purified  of 
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TABLE  IV.  This  table  provides  detailed  information  for  Fig.  2.  With  efficiency  r/=90%,  targeting  fidelity  F= 90%,  and  fiber  attenuation 
length  LaU= 20  km,  the  minimized  average  time  tavg  and  optimized  control  parameters  (L0,pc)  are  listed  with  respect  to  various  final 
distances  for  the  DLCZ  protocol. 


L  (km) 

20 

40 

80 

160 

320 

640 

1280 

2560 

5120 

10240 

tavg  (s) 

0.0013 

0.0065 

0.046 

0.44 

5.8 

95 

1900 

5.7  X  104 

3.0  X  106 

3.0  X108 

L0  (km) 

10 

20 

40 

40 

40 

80 

80 

80 

80 

80 

Pc 

0.26 

0.18 

0.13 

0.040 

0.011 

0.010 

2.7  X  10“3 

7.1  X  10“4 

1.8X  10“4 

4.6X  10“5 

Ffin 

90% 

90% 

90% 

90% 

90% 

90% 

90% 

90% 

90% 

90% 

vacuum 

and  multiexcitation  type  errors 

during  ENC.  The 

<<s2>  = 

2  DL0, 

(9) 

closest  analog  to  the  DLCZ  protocol  is  our  new  scheme 

without  ENP,  i.e.,  only  ENC.  At  longer  distances,  our  ap-  where  D  is  the  Phase  diffusion  coefficient  of  the  fiber.  If  the 
proach  is  further  improved  in  comparison  to  the  DLCZ  pro-  Phase  fluctuation  satisfies  Gaussian  distribution,  the  phase 
tocol  due  to  the  reduced  amplitude  of  vacuum  terms.  error  probability 


B.  Logical  errors 


P phase-err  \  ^ill  )  (1  e  °) . 


(io) 


So  far,  we  have  only  considered  the  effects  of  inefficiency 
that  maps  states  between  logical  and  nonlogical  subspace  by 
changing  the  number  of  excitations.  Besides  inefficiency, 
there  are  other  imperfections  which  preserve  the  number  of 
excitations  but  induce  errors  within  the  logical  subspace, 
such  as  interferometric  path  length  fluctuation  and  linear  op¬ 
tical  misalignment.  For  example,  the  interferometric  path 
length  fluctuation  leads  to  a  phase  difference  S,  which 
changes  |  '&ENG)a,b  into  a  mixture: 


For  example,  D=  10  3  rad2 /km,  L0=  1 0  km,  and  ppbase-err 
=  0.5%. 

Also,  a  small  probability  ( perr )  of  linear  optical  misalign¬ 
ment  per  ENC  or  ENP  step  is  modeled  as  depolarizing  errors. 
Later,  we  will  demonstrate  that  errors  within  the  logical  sub¬ 
space  restrict  the  final  fidelity  of  the  DLCZ  protocol,  while 
for  our  new  approach  additional  active  purification  can  cor¬ 
rect  such  logical  errors  to  achieve  high  fidelity. 


I VENG)a,b  -  ^(\H)a\V)b  +  eiS\V)a\H)b)  +  |tfV»ac)fc 

+  e2/<^^|vac)J//V>/,,  (8) 

where  (f>  is  the  static  phase  difference  between  the  left  and 
right  channels.  Since  the  last  two  terms  with  | HV)  will  be 
removed  during  the  first  level  of  ENC,  the  static  phase  4>  has 
no  effect.  However,  the  probability  of  being  in  an  undesired 
logical  state  j'L-)  is  sin2  |.  The  first  level  of  ENC  with  the 
combined  two  inputs  of  15y'_  and  fl,+  gives  <bA  producing  a 
phase  error  with  probability  2pphase.err=2(sin2  |),  propor¬ 
tional  to  the  variance  of  the  interferometric  phase  fluctuation. 
This  error  will  be  amplified  during  subsequent  ENCs  be¬ 
cause  the  survival  probability  of  the  state  <J?_  (the  logical 
error)  is  twice  as  much  as  that  of  <I>+  (the  desired  compo¬ 
nent). 

We  expect  that  there  is  little  correlation  in  phase  fluctua¬ 
tion  between  different  sections  of  the  fiber,  and  the  variance 
of  the  phase  fluctuation  is  proportional  to  the  length  of  the 
fiber 


IV.  SCALING  AND  TIME  OVERHEAD 
FOR  QUANTUM  REPEATER 

A.  Scaling  analysis 

Based  on  the  calculation  of  the  success  probability  at  each 
level  of  connection-purification,  we  can  obtain  the  estimated 
average  time  for  various  schemes.  In  Fig.  2,  we  compare  our 
approach  to  the  DLCZ  protocol  (see  Tables  IV  and  V).  For 
the  DLCZ  protocol,  the  average  creation  time  for  a  distant 
pair  contains  a  superpolynomial  contribution  (but  still  subex- 
ponentially)  with  distance,  due  to  instability  of  the  vacuum 
component.  For  our  new  scheme,  the  scaling  is  strictly  poly¬ 
nomial  with  distance  tavg^La,  where  the  exponent  a=a(  rj) 
explicitly  depends  on  the  efficiency. 

We  remark  that  the  DLCZ  protocol  is  slightly  more  effi¬ 
cient  for  short  final  distances  (L^  160  km).  The  DLCZ  pro¬ 
tocol  can  skip  the  entanglement  connection  and  exploits 
postselection  to  create  the  polarization  entangled  state,  while 
the  new  scheme  requires  the  first  level  of  ENC  to  eliminate 


TABLE  V.  This  table  provides  detailed  information  for  Fig.  2.  With  efficiency  y=90%,  targeting  fidelity  F=  90%,  and  fiber  attenuation 
length  Latl= 20  km,  the  minimized  average  time  tavg  and  optimized  control  parameters  ( L0,pc )  are  listed  with  respect  to  various  final 
distances  for  the  new  scheme. 


L  (km) 

20 

40 

80 

160 

320 

Lug  (s) 

0.0051 

0.020 

0.10 

0.68 

5.4 

L0  (km) 

5 

10 

20 

40 

40 

Pc 

0.26 

0.17 

0.11 

0.087 

0.037 

Ffin 

90% 

90% 

90% 

90% 

90% 

640 

1280 

2560 

5120 

10240 

45 

380 

3.3  X  103 

2.9  X  104 

2.6  X  105 

40 

40 

40 

40 

40 

0.017 

8.  IX  10“3 

4.0  X  10“3 

2.0X  10“3 

9.7  X  10- 

90% 

90% 

90% 

90% 

90% 
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FIG.  3.  (Color  online)  Average  time  versus  final  fidelity  with  no 
phase  errors.  For  fixed  final  distance  L=1280,  we  plot  the  (final) 
fidelity  dependence  of  the  (optimized)  average  pair  creation  time, 
for  the  DLCZ  protocol  (black  dashed  lines),  new  scheme  without 
ENP  (blue  solid  lines),  and  new  scheme  with  ENP  (red  dashdotted 
lines),  with  efficiency  77  to  be  90%  (circles),  and  95%  (no  circles). 

those  unwanted  components  and  prepare  the  polarization  en¬ 
tangled  state.  The  postselection  has  success  probability  1/2, 
two  times  more  efficient  than  the  success  probability  1/4  for 
the  first  level  of  ENC  of  the  new  scheme. 

For  the  new  scheme  (without  ENP),  we  can  use  the  stable 
probability  distribution  to  estimate  the  average  time: 

tavg  =  r0(L/L0)log2{1'5[2(2  -  <dV(3-2*»  (ll) 

where  h)~y~n~feL°ILatt  is  the  elementary  pair  generation  time, 
pc  is  the  elementary  pair  generation  probability,  L0  is  half  the 
distance  between  neighboring  repeater  stations,  and  L  is  the 
final  distance.  The  exponent  can  be  understood  as  an  over¬ 
head  from  the  finite  success  probability  for  ENC,  Penc 

rfil-lr/) 

~  2(2 -rjf  •  The  constant  1.5  in  Eq.  (11)  is  the  empirical  esti¬ 
mate  of  the  overhead  from  the  waiting  time  to  obtain  two 
independent  pairs  versus  the  single  pair.  In  Fig.  2,  the  differ¬ 
ences  between  the  simulated  data  and  Eq.  (11)  are  attributed 
to  the  empirical  factor  (1.5)  and  an  overall  factor  from  the 
overestimate  of  the  success  probability  for  the  first  level  of 
ENC.  According  to  Eq.  (7),  one  can  always  reach  good  final 
fidelity  if  the  elementary  pair  generation  probability  pc  scales 
as  L0/L.  Therefore  the  average  distant  pair  generation  time 

scales  exactly  polynomially  with  distance  tave^La,  with 
I  2(2-  Vf  \ 

«=l+log2(1.5)  +  log2(7(— j. 

B.  Comparison  between  different  schemes 

Besides  the  DLCZ  protocol  and  our  new  approach  with¬ 
out  ENP,  we  now  consider  a  scheme  with  ENP,  which  has 
one  phase-ENP  after  the  second  level  of  ENC.  We  may  com¬ 
pare  these  schemes  by  using  t-F  plots — a  parametric  plot  of 
tavg  and  F  as  a  function  of  excitation  probability  pc.  For 
given  noise  model  and  efficiency  77,  a  repeater  scheme  cor¬ 
responds  to  a  curve  on  the  t-F  plane. 

In  the  absence  of  interferometric  path  length  fluctuation 
(Fig.  3),  the  new  approach  without  ENP  is  about  five  times 


FIG.  4.  (Color  online)  Average  time  versus  final  fidelity  with 
phase  errors.  This  shows  the  same  optimized  approaches  as  Fig.  3 
but  with  finite  interferometric  path  length  fluctuation  characterized 
by  the  phase  diffusion  coefficient  D=10-3  rad2 /km.  Inclusion  of 
ENP  (red  dash-dotted  lines)  yields  a  dramatic  improvement  in  time 
for  high  fidelity  operations. 

faster  than  the  DLCZ  protocol,  for  77=  90%.  As  given  by  the 
previous  discussion,  this  improvement  is  due  to  better  con¬ 
trol  of  inefficiency-induced  imperfections.  There  is  a  time 
overhead  for  the  new  approach  (with  ENP)  as  compared  to 
the  new  approach  (without  ENP).  Within  each  implementa¬ 
tion,  the  higher  the  efficiency  77,  the  faster  the  quantum  re¬ 
peater.  For  high  final  fidelity  (l-T'^10%),  the  curves  ap¬ 
proach  straight  lines  with  slope  -1  because  t^p~l<x(  1 
-F)~\ 

When  interferometric  path  length  fluctuation  (leading  to 
initial  phase  error)  is  non-negligible,  active  ENP  is  needed. 
We  use  a  diffusion  model  for  the  path  length  fluctuation,  as 
detailed  in  Sec.  Ill  B.  In  Fig.  4,  t-F  curves  are  plotted,  as¬ 
suming  the  phase  diffusion  coefficient  D=  10-3  rad2/km,  cor¬ 
responding  to  pphase-err~  0.5%  over  L0=  10  km.  Unlike  Fig.  3 
where  only  inefficiency  is  considered,  there  is  an  upper 
bound  in  final  fidelity  for  each  implementation.  Both  the 
DLCZ  protocol  and  new  approach  (without  ENP)  suffer  from 
the  initial  phase  error,  with  final  fidelity  no  more  than  65%, 
while  the  new  approach  with  ENP  maintains  high  final  fidel¬ 
ity  up  to  97%.  For  high  retrieval  and  detection  efficiency 
( 77=  95  % ),  the  new  approach  (with  ENP)  can  produce 
1280  km  entangled  pairs  with  fidelity  90%  at  a  rate  of 
2  pairs/h,  even  in  the  presence  of  substantial  dynamical 
phase  errors. 

V.  OUTLOOK 

In  summary,  our  new  approach  to  long  distance  quantum 
communication  uses  a  different  qubit  basis  which  prevents 
the  growth  of  vacuum  and  multiexcitation  probabilities.  This 
keeps  the  ENC  success  probability  high  and  error  probability 
low,  and  leads  to  true  polynomial  scaling  even  in  the  pres¬ 
ence  of  realistic  inefficiencies.  We  can  achieve  a  bandwidth 
of  1  (or  2)  entangled  pair(s)  per  6  min  for  F~90%  (or  78%) 
and  negligible  initial  phase  error.  The  new  approach  also 
allows  active  entanglement  purification,  which  combined 
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with  built-in  purification  of  transmission  loss  errors  allows 
purification  of  arbitrary  errors  in  quantum  communication. 

Although  the  present  approach  shows  a  dramatic  im¬ 
provement  in  communication  rates  and  robustness  compared 
to  the  original  DLCZ  protocol,  the  bandwidth  remains  rela¬ 
tively  slow,  even  when  very  high  efficiencies  and  very  long- 
lived  quantum  memory  are  assumed.  While  such  high  effi¬ 
ciencies  might  ultimately  be  achievable  (see  Ref.  [23]  for 
recent  progress),  other  approaches  need  to  be  considered  that 
can  further  improve  the  effective  communication  bandwidth. 
For  example,  we  can  use  many  cells  per  node  to  improve  the 
bandwidth.  In  this  case,  the  improvement  is  at  least  linear 
with  the  number  of  cells,  making  it  possible  to  realize  long 
distance  (1280  km)  entangled  state  generation  bandwidth  of 
the  order  of  one  pair  per  second. 

A  simple  Monte  Carlo  optimization  of  efficient  use  of  the 
cells  shows  that  we  can  increase  the  bandwidth  by  a  factor  of 
r=M 112,  where  M  is  the  increment  of  the  factor  of  physical 
resources.  Recently,  it  has  also  been  suggested  that  multiple 
cells  can  be  used  to  further  facilliate  the  quantum  repeater,  in 
the  presence  of  memory  errors  from  the  quantum  memory 
[24,25], 

Note  added.  Since  completion  of  this  work,  a  work  by 
Zhao  el  al.  [26]  describing  a  similar  approach  to  improve  the 
DLCZ  protocol  has  appeared. 
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APPENDIX 

We  now  justify  two  main  claims  used  in  the  previous 
discussion.  (1)  For  the  DLCZ  protocol,  the  probability  ratios 
(m)/  (m)  ^  (m)  i  (m)  jncrease  wjt]1  the  nesting  level  m. 

(2)  For  the  new  scheme  (without  ENP),  the  probability  ratios 

PimJ p\'o  \c  and  P^mliti  P/'n^ic  remain  almost  independent  of  the 
nesting  level  m.  The  states  of  logical,  vacuum,  and  multiex¬ 
citation  types  depend  on  the  repeater  protocol.  For  example, 
the  entangled  logical  state  for  the  DLCZ  protocol  has  a  total 
of  one  excitation  stored  in  two  remotely  entangled  cells  [e.g., 
Eq.  (1)],  while  the  entangled  logical  state  for  the  new  scheme 
has  two  excitations  stored  in  four  cells,  with  one  excitation 
and  two  cells  on  each  side  [e.g.,  Eq.  (2)].  Thus  the  definitions 
of  p\'"\  ,  p'"],  and  are  different  for  the  two  schemes. 

1  logic ’  *  vac ’  1  multi 

APPENDIX  A:  NONLOGICAL  STATES 
FOR  DLCZ  PROTOCOL 

We  start  with  the  DLCZ  protocol.  First,  we  decompose 
the  density  matrix  (for  a  pair  of  distant  atomic  cells  x  and  y) 
into  components  with  different  excitation  patterns,  neglect¬ 
ing  the  interpattern  coherence 

Px,y  =  P oo77oo  +  Pio7rio  +  Pii7rn  +  Pioho  +  ‘  »  (Al) 

where  PijTTy  is  the  projected  density  matrix  to  the  subspace 
spanned  by  the  Fock  states  of  the  cells  {|f)A|y)v,|y)A-| with 


probability  p,j  and  normalized  density  matrix  7r,;-.  For  the 
DLCZ  protocol,  the  vacuum  type  of  state  is  7r00,  the  logical 
type  of  state  is  7r10,  and  the  rest  belong  to  the  multiexcitation 
type.  We  may  also  introduce  the  notation  corresponding  to 
Eq.  (5), 

Px,y  ~  PlogicPlogic  Pvac^uac  Pmulti77 multi  (  A2) 

with  pvac=p00 ,  Piogic=P to.  and  pmuiti=Pn+p2 o+‘  * '• 

Since  we  are  only  interested  in  the  coherence  properties 
for  the  logical  type  7 r10,  we  only  keep  track  of  the  probabili¬ 
ties  for  the  vacuum  and  multiexcitation  types  and  neglect 
their  coherences.  From  the  symmetry  for  the  two  cells,  we 
have 


iTji  =  (I0.v<*1  ®  \j)y(j\  +  \  j)x(j\  ®  \i)y(i\)  ■  (A3) 


In  the  rest  of  this  appendix,  we  will  add  a  superscript  to 
only  when  we  want  to  keep  track  of  the  coherence  for  that 
specific  term.  For  example,  indicates  that  there  is  coher¬ 
ence  between  the  states  |l)r[0)v  and  |0)v|  1)-^  after  the  mth 
level  of  ENC. 

Since  the  DLCZ  protocol  requires  that  the  probability  for 
the  multiexcitation  states  should  always  be  much  smaller 
than  the  probability  for  the  logical  states  (otherwise,  a  large 
fraction  of  multiexcitation  states,  accompanied  by  photon 
loss,  can  induce  significant  logical  errors),  we  regard  p\'"\ 
Pw(*Pw)  as  perturbations. 

We  denote  the  entangled  state  after  the  mth  ENC  as 


(m)  (m)  ,  (m)  (m)  ,  (m)  ,  (m) 

Px,y  =Poo  rroo  +  Pio  '7rio  +Pu  77  u  +  P20  ^20+  ■"  ■ 


(A4) 


We  now  connect  two  such  entangled  states  and 
for  the  (m+  l)th  ENC  via  the  superoperator  £ENC, 


(m 

xR,y 


A»'+ 1) 

lxr,-yR 


ZencIPx 


(A5) 


where  p"!+1)  is  the  unnormalized  density  matrix  for  the  en- 

xL’yR 

tangled  state  after  the  (m+  l)th  ENC.  Since  £ENC  is  a  linear 
operator  and  two  inputs  have  the  same  state. 


-'ENC 


P  n  T(  .  ^2  p pTT p 
a  p 


2  PaPp^ENC^-77^  77 p\  > 
a.p 


(A6) 


where  ^c[ tra,  1 rj  =  \£ENdTTa,  1 r^]  +  ^enJ-^P’  77 J.  Now 
we  calculate 


for 


For  example, 


£mct7ra’  77 p\ 

{(mi)  1 

7700>7710  ’7711’7720>  •••  )• 


r  (rn)  (m)-l  _  V  (rn+1)/  ,  Ilil - ^  (A7a) 

L^io  ’^to  J  _  ^to  +  9 


csym 

C-ENC  L  "  10  >  "10  J  _  "10 


CT77ioUoo]  =  f77oo,  (A7b) 
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^cKioUi]  =  “''Hi  +  ’X1  -  ^io. 


r  M  ^-rL  , 
ifiJVcL^lO  ,W2oJ-  4^20+  2  77-10 


3  77(1  -  77) 2 


■^oo. 


fsvc[7roo>'^oo]  =  0. 


the  multiexcitation  state  (iru)  are  at  least  doubled,  relative  to 
(A7c)  the  logical  state  7rj™+l\  after  each  ENC. 

In  terms  of  the  notation  corresponding  to  Eq.  (5)  (pvac 
=Poo ,  Piogic=P  10.  and  pmm=pn+p2 0+*").  we  have 


(A7d) 

(m+l)  ~(m+ 1)  On) 

P  vac  P  00  ^  ^  7  vac 

(m+l)  —  ~(m+ 1)  ^  (m)  ’ 

P  logic  P 10  P  logic 

(Alla) 

(A7e) 

~(m+ 1)  (m) 

P multi  _  ^  P multi 
~ (m+l)  ^  ( m)  ' 

Plogic  Plogic 

(Allb) 

^HVc['7rOO>77'll]  =  (A7f) 

£Senc\-  ^  ^o]  =77(1  —  77)  7700 .  ( A7  g) 

The  logical  type  of  state  after  the  (7«+l)th  ENC,  is 

the  average  of  the  states  and  7r10,  with  relative 

weights  V10 Vio’  and  2 77Q  -  v)PwPn)+  vk 1  “  v)PmPw 
+  2  r/PooPn  \  respectively. 

Then  we  calculate  the  (unnormalized)  density  matrix  after 
the  (m+l)th  ENC, 


The  ratio  of  nonlogical  states  to  logical  states  is  at  least  dou¬ 
bling  with  distance.  As  discussed  in  the  main  text,  it  is  these 
unstable  nonlogical  states  that  leads  to  the  superpolynomial 
scaling  for  the  DLCZ  protocol. 

APPENDIX  B:  NONLOGICAL  STATES 
FOR  NEW  SCHEME 

For  the  new  scheme,  we  can  similiarly  decompose  the 
density  matrix  [following  Eq.  (5)] 


=  MoWo  +  Mo+1Mo+I)  +  Mr  %  1 

+  Mo+1)'7720  + 

with  the  probability  coefficients 


(A8) 


Am+l)  __  (rn)|  (m)  (1 _ V}_  (m) 

P  00  VP  10  1 P  00  +  P 10 


(A9a) 


Am+l)  __  Hn(m)  (m) 
P 10  ~  2  ”10  "10  > 


P 11  vPioPn 


Am+i)  V  (m)  (m+1) 

P 20  ^  PlO  P20 


(A9b) 

(A9c) 

(A9d) 


P  -P0077b0  +  Pl0'7r10  +  Pll'7rll  +  P  10^10  +  Pll'Irl\  +  P307r30"  ’  . 

(Bl) 

where  /7(-;7r,-;  is  the  projected  density  matrix  to  the  subspace 
spanned  by  the  Fock  states  with  i  (or  j)  photons  in  a  cell-pair 
and  j  (or  i)  photons  in  b  cell-pair,  with  probability  ptj  and 
normalized  density  matrix  77 y,-.  For  the  new  scheme,  the 
vacuum  type  of  states  consists  of  7 r00  and  7r10,  the  logical 
type  of  state  is  7rn,  and  the  rest  belong  to  the  multiexcitation 
type. 

We  use  a  perturbative  approach,  by  assuming  Poo7Tbo> 
/j k)7T|0,  and  p\]TT\  \  are  the  dominant  terms,  and  the  rest  terms 
are  perturbations  of  order  pc  (terms  not  listed  are  of  order 
pi).  We  eliminate  those  irrelevant  perturbation  terms  (e.g., 
p307T30)  because  after  one  level  of  ENC,  they  are  suppressed 
to  0(pl). 

Suppose  the  entangled  state  after  the  /nth  ENC  is 


to  the  leading  order  with  respect  to  the  perturbations  of  M'!’* 
and  p2Q’. 

Finally,  we  divide  these  probabilities  by  /L™+11,  and  obtain 


McT"  JpS  (1-17)1.  .fg 
Mr"''  ^  )  psr 


(A10a) 


— Tm+l) 

P 11 


PlO 


Jm) 

??1L 

~  On)  ’ 
PlO 


(A  10b) 


~(m+l) 

PlO 

={m+ 1) 
PlO 


(m) 

P20 

„(m)  ' 
PlO 


(AlOc) 


Since  the  normalization  does  not  change  the  relative  ratio 
between  the  probabilities,  the  above  perturbative  estimate 
tells  us  that  the  fractions  for  both  the  vacuum  state  (7r00)  and 


P$  =  pftVoo  +  pftV  10  +  M7M'  f  +  pioV20.ll  +  pi ol  ^20, x 
+  P2U7T2i,ll  +  piT.x^i.i.  +  '  •  ‘  ■  (B2) 

Notice  that  we  need  to  distinguish  two  possible  types  of 
states  for  tt2q  because  they  behave  differently  during  ENC. 
The  first  type  of  states  (denoted  as  7720,11)  has  both  photons 
stored  in  the  same  cell,  and  after  retrieval  the  photons  will 
have  the  same  polarization,  follow  the  same  path  way,  and 
trigger  the  photon  detector(s)  on  the  same  side  of  the  PBS 
[Fig.  1(b)].  The  second  type  of  state  (denoted  as  7720,1)  has 
two  photons  stored  in  different  cells,  and  after  the  retrieval 
the  photons  will  have  orthogonal  polarization,  split  at  the 
PBS,  and  trigger  photon  detectors  on  both  sides  of  the  PBS 
[Fig.  1(b)].  Thus  the  second  type  of  state  is  more  likely  to 
give  the  correct  click  pattern  and  thus  propagate  the  error  to 
the  next  level  of  ENC.  Similarly,  we  introduce  7r21  h 
and  7r2iix- 
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In  the  rest  of  the  discussion,  we  still  follow  the  conven¬ 
tion  from  the  previous  appendix  that  77, ;  =  tt;,  and  we  will  add 
a  superscript  m  to  irlJ  only  when  we  want  to  keep  track  of  the 
coherence  for  that  specific  term. 

We  now  connect  two  such  entangled  states  p"l\  and 
pm\  for  the  (;n+l)th  ENC, 


— {m+ 1)  c  V  {m)  ( m )  i 

Pa,,bD  ~  ^ENG-Pa,  ,b^  Par,,bB\'> 


(B3) 


where  ’p'"+h "  is  the  unnormalized  density  matrix  for  the  en¬ 
tangled  state  after  the  (in  +  1  )th  ENC. 

Now  we  calculate 


^'EVc[7ra’7rJ s]  —  2^'£A'c['77a>7r/3]  +  l^ENci77/}’  77  J 


for 


77a>77/8  £  Wo0>7710’77n*’7720,ll(or  _L)^ '7T'21,|l(or  X)>  ■■■  }• 

For  example, 

9 

csytn  r  (m)  (m  h  7 T  (m+ 1 )  r 

-WC77  n  .77ii  ]  =  y77ii  , 


csvm  r  (m)  *i  'i 

‘-'HVCX77!!  ’  '7r toJ  —  ,  ^lO, 


jt. 

4 


£S1ncW nUoo]  =  0, 

=r(  1  -  pHio, 


nsym  r  (m)  n  _  rf  ,  rf(*  ~  rf 

^ENO-^n  ’7721,<5.J  ~  ^^21,3+  0 


s  m  rf 

ZenJ-^W’  77  to]  =  -g-77!®’ 

^enJ-17  io>  ^oo]  =  0> 


(B4a) 

(B4b) 

(B4c) 

(B4d) 


^hvc[7Tio>77'21,ii]  : 


rf(  1  -  V)  rf 


£svc[7710>7721,x]  -  ^771,1  + 


V 2  rf(^~v)  rf 


^vc^oo^ixhO, 


“^EVcI77 00’  7r21,x]  —  ,,  TUo- 


_7710  +  o  ^O.X* 

(B4i) 
(B4j) 

(B4k) 
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The  logical  type  of  state  after  the  (m+l)th  ENC,  rf™+1\  is 
the  average  of  the  states  7r<1"1'+1^  and  7ru,  with  relative 
weights  fp'/'/Vi'"’  and 
respectively. 

Then  we  calculate  the  (unnormalized)  density  matrix  after 
the  (m+  l)th  ENC, 

~(m+l)  _  ~(m+ 1)  ~4m+ 1)  ,  ={m+l)  (m+1)  ,  -Am+1) 

Pa,.bK-P  00  ^OO  +  PlO  ’Uo  +  Pll  ^11  +  P20.II  ^20,  II 


,  (m+1)  ,  (m+1)  ,  (m+1)  , 

+  P20.X  ^20, X  +  P21JI  7721,|I+P21,X  ’Dl.X  + 


with  the  probability  coefficients 

~(m+l)  _  7  (m)  (m) 

Poo  TT^io  Pio  > 


2\+2(l-^1T)(pW+^) 


~(m+l)  HLn(m)  (m) 
Pl0,S  ^  PlO  P2i,a> 


fr(m+l)  V_  (m)  (in) 

Pll,S  Pit  Pll,S- 


(B5) 

(B6a) 

(B6b) 

(B6c) 

(B6d) 

(B6e) 


Finally,  we  divide  these  probabilities  by  p"rf\  and  to  the 
leading  order  (i.e.,  the  zeroth  order  of  pc)  we  obtain 


(B4e) 

pti+l)  ~  i/p^v 

MT+1)  4  \pfrV  ’ 

(B7) 

(B4f) 

Mo+1)  _  p^ 

fjfm+l)  „(m)  ’ 

Pll  Pit 

(B8) 

(B4g) 

~(m+l)  (m)  (m) 

^20,5  LP 10  P2\,S 

~(m+l)  9  _(m)  (in)  ’ 

Pll  ZPll  Pll 

(B9) 

(B4h) 

~<m+l)  (m) 

Pll, 8  _  ^21, <5 

— (hi+1)  (m)  ' 

Pll  Pll 

(B10) 

Furthermore,  all  these  ratios  remain  constant  (to  order  pc), 
which  justifies  the  claim  that  the  probabilities  for  different 
types  of  states  remain  stable  for  all  higher  levels  of  ENC. 

If  we  further  introduce  rf"1* 1  f  =  pZ\+ 1  *  +p!rf '  ^  p\'"+1) 

— (m+1)  ,  ~<m+l)  v  Am+irUm+l)°°  /10 

=  P It  •  and  Kiulti  =^S=l±P20,S  +Pl\,S  >  We  haVe 


logic 


fplm>^[l  +pt,nl  +  0(pc)l  (B 1 1  a) 


p£l)~Jp£luP%&l  +  0(pc)], 


(B  lib) 
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iffi?  ~  y +  0(Pc)l  (Bile) 

where  pim+1)  is  the  probability  for  the  new  logical  error  from 

1  err,new  sr  j  a 

the  multiexcitation  states  (accompanied  by  photon  loss) 


(,„+p  „ 

P  err.new  (m)  (m) 

P 11  P 11 

(B12) 
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